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YONSEI EVOLUTIONARY POPULATION SYNTHESIS (YEPS) MODEL. I. 
SPECTROSCOPIC EVOLUTION OF SIMPLE STELLAR POPULATIONS 

Chul Chung Suk-Jin Yoon^ '^, Sang-Yoon Lee^, & Young- Wook Lee^ 

ABSTRACT 

We present a series of papers on the year-2012 version of Yonsei Evolutionary Population 
Synthesis (YEPS) model which is constructed on over 20 years of heritage. This first paper 
delineates the spectroscopic aspect of integrated light from stellar populations older than 1 Gyr. 
The standard YEPS is based on the most up-to-date Yonsei-Yale stellar evolutionary tracks and 
BaSel 3.1 flux libraries, and provides absorption line indices of the Lick/IDS system and high- 
order Balmer lines for simple stellar populations as functions of stellar parameters, such as 
metallicity, age and a-element mixture. Special care has been taken to incorporate systematic 
contribution from horizontal-branch stars which alters the temperature-sensitive Balmer lines 
significantly, resulting in up to 5 Gyr difference in age estimation of old, metal-poor stellar 
populations. We also find that the horizontal branches exert an appreciable effect not only on 
the Balmer lines but also on the metallicity-sensitive lines including the magnesium index. This 
is critical to explain the intriguing bimodality found in index distributions of globular clusters 
in massive galaxies and to derive spectroscopic metallicities accurately from various indices. A 
full set of the spectroscopic and photometric YEPS model data of the entire parameter space is 
currently downloadable at http://web.yonsei.ac.kr/cosmic/data/YEPS.htm, 

Subject headings: stars: general — stars: abundance, evolution, horizontal-branch — globular 
clusters: general 

1. INTRODUCTION 

The Evolutionary Population Synthesis (EPS) technique is a key tool for interpretation of integrated 
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2011;Pforr et al.ll2012h . Combined with recent development in high precision observations, the EPS models 



are becoming more important for the analyses of various stellar populations in galaxies. 

In a series of papers, we intend to present the Yonsei Evolutionary Population Synthesis (YEPS) model 
for the spectroscopic and photometric evolution of simple stellar populations (SSPs). This paper, as the 
first paper of the series, describes the spectroscopic aspect of our YEPS model. The model is constructed 
by the YEPS Fortran code package, which has been improved and exploited for the p ast 20 years by 
many studies related to (a) synthetic color-magn i tude diagrams for individual stars (e.g.. 
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2012). The forthcomina Paper II 



(Yoon et al. 2012, in preparation) will present the photometric evolution of stellar populations. The latter 
papers in the series will discuss the effect of the different choice of model ingredients and input parameters 
on the model, as well as the application of the YEPS for the early-type galaxies that have composite stellar 
populations. 

The standard YEPS model has been constructed based on the Yonsei- Yale (Y^) stellar evolution models 
(Kim et al. 2002; S. Lee et al. 2012, in preparation) and th e BaSeL flux library (Westera et al. 2002). The 
absorption-line index model employs the Lick/IDS system (IBurstein et al.ll 19841 ; iFaber et al.ll 19851 ; IWorthey 



19941 ; IWorthey & Qttavianilll997l ; ISchiavonll2007r) . which defines 25 absorption lines produced by various 
elements at the surface of stellar atmosphere. The Lick/IDS system u ses the spectra of nearly 460 stars to 
cover a wide parameter space of t emperature, gra vity, and metallicity feuzzoni et al. 1992 , 1994 ; Worthey 



19941 ; IWorthey & Qttavianil 1 19971 ; ISchiavonI 120071) . The system, however, does not consider the grid of a- 



elements enhancem ent, and thus the en hancement should be treated theoretically. We a pplied the g- e lemen t 



correction ter r ns by iKom 



etal 



Thomas et all (l2003h . and lSchiavonI (120071') . 



(I2OO5I') to our model, following the schemes used by 



Trager et all (l200nh . 



The YEPS model has been built with particular interest in dealing with the core helium burning 
horizontal-branch (HB) stars. Since the presence of hot stars in globular clusters (GCs) and galaxies changes 
the overall shape and absorption feature of their spectral energy distributions (SEDs), especially at short 
wavelengths, the impa ct of hot HB stars (> 8000 K) has been a topic of great interest in the EPS community 



over the p ast 20 ye ars (Lee et al.^l 199d. I^1994J ; IWorthe 



1994 



3Em 



Lee et al 



200c 



; iThomasj 



etal.ll2003l ; lLee etal. 



Lee et all (|2000) first demonstrated that H/3 



|2005b; S chiavonl l2007; Yoon et al. 2006, 2008L l2011a 
absorption index — th e most popular age indica tor — is significantly enhanced by the presence of blue HB 
stars. More recently, Yoon et al. (2006, 2011al jb[) show that the systematic metallicity-dependent variation 
in HB temperature leads to the nonlinear relationship between metallicity and broadband optical colors. 
Despite the fact that hot, blue HB stars exert a strong effect on properties of integrated light from GCs and 
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galaxies, most EPS models to date took into account the details of HBs in a fairly limited manner. The 
YEPS, by contrast, elaborates the HB effect not as a merely contamination factor but as a crucial part of the 
EPS model for various spectroscopic and photometric observables. 

The paper is organized as follows. Section 2 describes the constructing procedure of the YEPS model. 
Section 3 presents the results of our stellar population simulations and the comparison of our model with 
observations. Section 4 discusses the implications, and finally Section 5 summarizes our results. A full 
set of the spectroscopic and photometric YEPS model data of the entire parameter space is available at 
http://web.yonsei.ac.kr/cosmic/dataA^EPS.htm 



2. CONSTRUCTION OF THE YEPS MODEL 

The YEPS model provides, for a given stellar system, (a) the synthetic color-magnitude diagrams 
(Section 2. 1), (b) the synthetic spectral energy distributions (Section 2.2), (c) the integrated absorption line 
indices (Section 2.3), (d) the integrated magnitudes and broadband colors, and (e) the integrated surface 
bright fluctuations. Tabled] summarizes the ingredients and input parameters of the YEPS model. 



2.1. Synthetic Color-Magnitude Diagrams 



The standard YEPS model is constructed based on the most up-to-date Yonsei-Yale (Y^) stellar evolu- 
tionary tracks. For the evolution ary phases from the main sequence (MS) to the tip of the red giant branch 
(RGB), we used Y^-isochrones (IKim et al.ll2002l : Lee et al. 2012, in preparation), covering the metallic- 
ity grids from Z = 0.00001 to 0.08 with three different values for the a-elements enhancement ([a/Fe] 
= 0.0, 0.3, and 0.6). Th e mixture pattern of a-elements enhancement in Y^-isochrones follows that of 
VandenBerg et al.l (120001 ). The Y^ stellar evolutionary libraries adopt the galactic helium enrichment pa- 
rameter of AF/AZ = 2.0 with the primordial helium abundance of F = 0.23. To examine the effect of the 
different choice of the evolutionary tracks, we comparatively used the BaSTl stellar evolutionary tracks 
(IPietrinfemi et al.ll2004l) with metallicities from Z = 0.0001 to 0.04 for the two a-elements enhancement 
cases ([a/Fe] = 0.0 and 0.4). The Y^ stellar libraries include helium diffusion and BaSTl stellar libraries 
include the atomic diffusion of both helium and metals. As will be demonstrated below, the major features 
of our model do not depend on the specific choice of stellar libraries. 

We adopt the Salpeter initial mas s function (IMF ) for our standard set of simulations to assign the 
number of stars along given isochrones. IWortheyl (119941) presented the generahzed Salpeter IMF of the form 

dN _ Mtot{\-x) 
dM~M}. 



M" 



-(1+^) 



(1) 



where dN is the number of stars within the fixed mass bin dM, and M/ and M„ are the lower and upper 
mass cuts, respectively. From this, we calculated the IMF of an SSP that consists of an single-metallicity 
and single-age population. We have applied 10^ stars within whole mass range of IMF (from 0.2 to 5.0 
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Mq). We adopted the standard Salpeter index (x = 1.35) over the whole mass range. The choice of x exerts 
a fairly small effect on the overall shape of UV-to-IR SEDs and hence on broad-band colors and absorption 
indices (IPark & Led 1 19971) . This is more so for old stellar populations for which massive stars are already 
evolved off the MS because the index x controls the fractional contribution from the massive stars. H owever, 
it is n oteworthy that x leads to significant variations in total absolute magnitude of the model SSPs (ITinsley 
19721 ) because the absolute flux level of SEDs is a function of the total stellar mass. 

For the synthetic HB modeling, we used Y^-HB tracks (Lee et al. 2012, in preparation) that are fully 
consistent with Y^-isochrones in terms of the input physics and assumed parameters. The Y^ HB tracks 
cover a wide range of HB total mass from 0.4438 M© for Z=0.06 and 0.5037 M© for Z=0.0001 to 1.5 Mq 
for all metallicities to incorporate the wide variation of HB morphology. In order to simulate the mass 
dispersion of HB stars, we have used the Gaussian HB mass distribution of the form 



P(M) (X exp 



-{M-{Mhb)? 
2al 



(2) 



where P(M) is the probability density function of the HB mass, and (Mhb) is the mean m ass of the HB at 
a give n metallicity and age. The standard model assumed the value of to be O.O15M0 (|Lee et al.lll990 . 
19941) . On average, the number of HB stars at given metallicity and age is 350 in a single simulation. In 
addition, in order to avoid small number statistics in HB modeling, we have repeated the simulation 10 times 
to get the averaged continuum flux at given metallicity and age. 

Figure[T]shows how the HB morphology of YEPS modeflis calibrated to the observations. The HB type 
is defined as (B-R)/(B + V +R), where B, R and V are the numbers of blue and red HB stars and RR Lyrae 
variable stars, respectively (|Lee et al.lll994|) . Filled circles represent the oldest inner-halo (Rcc < 8 kpc) the 
Milky Way globular clusters (MWGCs), and open circles (8 < Rqc < 40 kpc) and triangles (Rcc > 40 kpc) 
represent the outer-halo MWGCs. Solid lines from top to bottom are the HB type variation of the YEPS 
model with varying ages. The free parameter, the Reimers r nass-loss efficie ncy parameter r/, is used to 
calibrate our model HB t ypes to the observations. W e adopted iReimersI (|l977l )'s empirical formula for the 



mass loss along the RGB ( Rood|[l973 : Lee et al.lll99C ). The formula takes the form of 



■ dM 



dt ocr?^, where L,g, 

and R are luminosity, gravity, and radius of stars, respectively. The comparison of models and observations 



suggest a 1] of . 63 under the assumpt i on that the mean age of i nner-halo GCs is 12 Gyr (IGratton et al. 



19971 : lReidlll997l : IChabover et al.lll998 



Marm-Franch et al 



20091 : 



Potter et al 



20101). Figure |2] shows the 



color-magnitude diagrams (CMDs) for selected model SSPs (red and blue arrow s in Figure H]). I n general, 
the HB type becomes redder with increasing metallicity and decreasing age (e.g. jLee et al.ll 19941 ). 



'Recent observations and modeling indicate that tlie abundance anomaly, especially in He and CNONa, is present in the MWGCs 
with multiple stellar populations. Although the variation in He and CNONa abundance among GCs in the MW is large, the average 
variation between GC systems in different galaxies, as a whole , is not expecte d to vary greatly. Additionally, only 30 % of the 
MWGCs are significantly affected by the enhancement in He ( Lee et al.ll2007n . and this suggests that HB morphologies in the 
majority of the MWGCs are mostly controlled by total metallicity and age. However, if the average He enhancement in the 
MWGCs is not archetypical, and it varies significantly from one galaxy to another, our models presented here would need further 
revisions to reflect this. 
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2.2. Synthetic Spectral Energy Distributions 



Spectral energy distributions (SEDs) of SSPs are generated based on the synthetic CMDs (Section 2.1). 
The CMDs give the stellar parameters of individual stars in given SSPs, including effective temperature (T), 
surface gravity (g), global metalUcity ([Z/H]), and a-element enhancement ([a/Fe]). To derive theoretica l 
spectral fluxes (in units of [erg/s/cm^/ AD, we use the spectral library of BaSel 3.1 (IWestera et all 120021) . 
BaSel 3.1 is based on the expertise of Kurucz ( 1992t) and BaSel 2.2 ( Lejeune et al. 1998 ). and provides 
extensive and homogeneous grids of theoretical flux distributions calibrated to the colors of the MWGCs 
at all levels of metallicity. The library covers effective temperatures from 2,000 K to 50,000 K, gravities 
in a solar unit from logg of -1.02 to 5.50, and metallicities [Z/H] from -2.0 to 0.5. Note that the BaSel 
3.1 library assumes scaled-solar a-elements. We thus choose to use the total metallicity [Z/H] for the 
construction of SEDs of SSPs, rather than the iron abundance [Fe/nil. The equation [Z/H] = [Fe/H] + 
A[a/Fe] relates [Z/H] to [Fe/H]. In our model for [a/Fe] = 0.3, the factor A which depends on the a- 
element mixture of the model equals to 0.723. 



We calculate the expected flux Fx at a distance d using the form 

L 1 



F\ = 47r X 



(3) 



where Hx, L, Tgff and a are the flux intensity, luminosity, effective temperature of a star, and the Stefan- 
Boltzmann constant, respectively. From these fluxes of individual stars, integrated fluxes of all stages of 
stellar populations in the synthetic CMDs — MS to RGB, HB and post-asymptotic giant branch (PAGB) — 
are calculated using the following summation form: 



r'total 



: pMS _^ pRGB _^ pHB _^ pPAGB 



(4) 



Total mass of SSPs at given age and metallicity is normalized to 10 Mq. 



2.3. Absorption-line Strength Indices 



The absorption-line indices of the YEPS model are calculated using the polynomial fitting fiinctions. 
The fitting functions are derived from spectra (in the 4,000-6,000 A region) of Galactic stars and yie ld the 
line strengths as functions of stellar atmospheric parameters — metallicity, temperature, and gravity (IRose 



1985'; 'Jones 
1992. 1994 



1995; 



Vazdekis 



1999; 



Fabe?'^ 



1973; 



Rose 



Worthey & Ottavianilll997 ; Schiavon 



200' 



19841 : iDiaz et al.lll989l : IWorthevlll994l : iBuzzoni et al 
. The empirical polynomial fitting functions, com- 



bined with continuum levels of model SEDs, generate the absorption-line indices of SSPs. 



^In some photometric broadband colors (e.g., U —B), broadband colors for g-enhanced mix ture are better reproduced by scaled- 
solar spectra with the same [Fe/H] of the a-enhanced mixture, not [Z/H] i Cassisi et al J2004 ). 



-6- 



For the standard YEPS absorption index model, we use IWortheyl (119941) and IWorthey & Qttaviani 



(Il997h 's (hereafter W94 and W97) polynor nial fitting functi ons for 25 absorption indices of the LicMDS 
system. As a compa rison rnodel, we adopt ^ Schiavon (2007)'s (hereafter S07) fitting functions in the blue 
wavelength based on I jones (1999) stellar library. Figure [3] shows an example of the fitting functions (H/3 
line) as given in W94 and S07. The two sets of fitting functions agree well with each other. We note that the 
H/3 fitting function of SO? exhibits greater meta llicity-sensitivity for giant stars than that of W94 by virtue 
of a more recent spectral library by ljonesi (|1999h . 



From these fitting functions, we derive the equivalent width of a single star at a given temperature, 
gravity, and metallicity. The empirical fitting function provides an absorption index (/), which is transformed 
to the equivalent width (EW) using the following equations; 



EWO = fcX 



1-7 
AA 



EW(mag) = fc X 10 



-0.4/ 



(5) 



(6) 



where AA and fc are the index bandpass and flux continuum of an SED for each index. The continuum 
fluxes (fc) of each index were taken from the definition of Lick standard system described in W94. The unit 
of mag is used for CNi, CN2, Mgi, Mg2, TiOi, and TiOi lines. 

After determining the EW and fc of individual model stars, we finally calculate the integrated indices 
by summing the continua and EWs of all stellar populations of SSPs using the following formulae. 



Integrated 
i'^ Index 



AA' X 



\ E;[/C..,X(1-^; 

T,jfc.j 



-2.5log 



12jfc.j 



(7) 



where fcj, Ij, and AA' are the / model star's continuum, absorption index, and the width of the i index 
bandpass, respectively. 



2.4. Treatment for the Enhancement of a-elements 



The absorption indices for the case of enha nced g-elements are modeled as follows. We use stellar 
evolutionary tracks with enhanced a-elements ( Kim et al.ll2002h . For the stellar atmosphere model, we 
adopt the a-elements mixture ratios that is identical to that of stellar tracks for the sake of consistency. 
Table |2] shows the a-elements mixture of the Y^ stellar evolutionary model and BaSTI, as compared with 
the scaled-solar abundance ratio of metals taken from iGrevesse & Noelsl (|l993h . Our YEPS models for 
enhanced a-elements follow the two a-element mixtures of Y'^ and BaSTI. 



We then applied the response functions of a-elements by iKom et all (120051) (hereafter K 05) to apply 
the a-elements mixture to YEPS. Compared to the previous work by iTripicco & Belli (| 19950 . K05 have a 
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more extended metallicity space ranging from [Fe/H] = -2.25 to +0.67, and provide the response functions 
for 25 Lick absorption indices for three evolutionary phases (dwarfs, turn-offs, and giants). The response 
functions that K05 provide are the first partial derivative dl/d{Xi\ of the Lick index /q when an abundance 
increment of two times logarithmic i'^ a-element (C, N, O, Mg, Fe, Ca, Na, Si, Cr and Ti) is assumed. As 
described in TMB03, it is appropriate to expect / oc exp([X,]) for the optimal approximation. Hence, the 
Taylor expansion for In / instead of / is an adequate approach for the variation of the Lick absorption indices 
due to a -element abiindance changes. Neglecting the higher-order derivatives and following the notation 



/?o.3(0 of iTrager et all (120001) . we can express the Taylor expansion in the following forms, 



where Ro.sii) = l//o x dI/d[Xi] x 0.3 is the K05 index response for increased a-element / by 0.3 dex, and / 
and /q are the absorption index before applying the K05 response function and the model absorption index 
of the K05 parameter space, respectively. The exponential scale of Eqn. [8] yields 

Inew=llle^piR0.3im °' ^- (9) 

(=1 

Fitting functions of W94 and S07 give negative values for absorption indices when stellar populations 
are young, metal-poor (e.g., CN, Ca, and Fe lines) or old, metal-rich (e.g., H/3, H7, and ¥16). Since our 
calculation of the a-element fractional change is based on the logarithmic Taylor series of In/, we must avoid 
negative values of the absorption indices. Table [3] lists the negative minimum values of YEPS absorption 
index models when we adopt W94 fitting functions together with and BaSTI libraries. The simplest 
way to correct negative values in the fractional index change is to shift the negative indices into the zero or 
positive value and then compute the fractional change. We used the correction term 6 as listed in Table [3l 
and applied the following equation to correct absorption indices with negative values. 



/.».-« = (/-«)n«p(;^a[|o3) ■ (10) 



After the derivation of 1-6, we scale back the index Inew-^ 



3. RESULT OF STELLAR POPULATION SIMULATIONS 

3.1. The Effect of HB Stars on Simple Stellar Populations 

Figure |4] shows the most temperature-sensitive indices (Balmer lines, CNi, and G4300) as a function 
of [Fe/H] and highlights the effect of HB stars on integrated absorption line strengths of SSPs. In the YEPS 
model without HBs, H/3 index decreases monotonically with increasing [Fe/H] at given age. This trend 
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is due solely to the temperature variation of turn-off (TO) and RGB stars, which becomes cooler as the 
metallicity of stellar population increases. By contrast, in the models with HB stars (solid lines), the H/3 
and H^f strengths are significantly enhanced by the presence of blue HB stars at the metal-poor regime. 
CNi and G4300 are also sensitive to hot HB stars (W94) in the sense that these lines get weakened by 
the presence of hot HB stars. The models with even hotter HB stars in the most metal-poor GCs tend to 
approach the models without HBs, as the hotter HB stars are dim and too hot to have significant effect on 
these indices. 

The fact that blue HB stars can mimic young, hot TOs in the integrated fluxes of GCs has significant 
implication for determining ages and metallicities of stellar populations. For instance, when the H/3 strength 
is used as an age indicator, a 7 Gyr model with no HB stars exhibits H/3 strength identical to a 12 Gyr model 
with HB stars at [Fe/H] ~ -1.6. This suggests that one could seriously underestimate age of stellar popu- 
lations for old (> 10 Gyr) populations. Also, HB stars affect other absorption indices known as metallicity 
indicators (e.g., Mgb, Fe4383, and (Fe)), although the amount of change is small compared to the Balmer 
indices (see Section 4 below). 

Given the wide implication of the HB effect, it is highly important whether or not such an effect depends 
on the specific choice of stellar libraries for evolutionary tracks and the fitting functions. In this study, we 
employ two stellar libraries (Y^ and BaSTI) and two empirical fitting functions (W94-I-W97 and S07). The 
four combinations of the models are presented in Figured Solid and dashed lines represent the H/3 strengths 
for the models with and without the HB prescription. Comparison shows that all models agree well for the 
H/3 index. For instance, all the 12 Gyr models with HB stars show the same amount of H/3 enhancement 
of about 0.6 A at [Fe/H] ^ -1.6 comp a red to the model without HB stars. Our additional test using the 



high-resolution spectra of iMunari et al.l (120051) also confirms that, compared to the model without HBs, the 
model with HBs at the same condition shows enhanced H/3 by about 0.6 A. We note that comparison between 
models without HB stars (dashed lines) shows that the stellar library produces a larger gap between iso- 
age lines than the BaSTI library. The reason for this is that the Y^-isochrones have larger TO temperature 
gaps than the BaSTI isochrones. In addition, the fitting functions of S07 yield slightly weaker H/3 indices in 
the metal-poor regime than W94. 



3.2. The Effect of a-elements on Simple Stellar Populations 

Figure ^demonstrates the effect of a-elements on Lick indices using two indices (H/3 and Mg b) that are 
sensitive to the effective temperature and a-elements, respectively. The middle column shows the effect of 
a-elements on H/3 as functions of [Fe/H] and age. The strength of H/3 index without HB stars (dashed lines) 
decreases with increasing [a/Fe]. This is because temperature of TO stars decreases with increasing [a/Fe]. 
Hence, the use of stellar libraries that incorporate enhanced a-elements is crucial to predict accurately the 
strength of H/3. The total Z increment due to the a-elements enhancement also exert effect on the H/3 
index with HB stars. The way the a-elements affect the HB types is shown in the left column of Figure |6] 
As [a/Fe] increases, the models with enhanced [a/Fe] (blue lines) show redder HB types compared to the 
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scaled solar models (green lines) at given [Fe/H]. The strengths of H/3 index in the middle column reflect 
the trends on the HB types shown in the left column. Interestingly, the 15 Gyr model (solid lines in the 
left column) has blue HB type at [Fe/H] = -2.5 ~ -1.5, but the contribution of blue HBs to H/3 is small 
compared to the 12 Gyr model at same metallicity. The first reason for this is that the strength of H/3 reaches 
its maximum at T^ff ~ 9,500 K. So, the contribution from hotter (T^ff > 9,500 K) HBs in older and/or more 
metal-poor populations to the K/3 absorption becomes smaller. The second reason is that those HBs with 
hot temperature are fainter. Hence, their contribution to luminosity-weighted absorption indices becomes 
relatively small. 

It is important to note that Mgb, the well-known tracer of [Fe/H] and [a/Fe], is also affected by the 
variation of HB morphologies in their absorption strengths (see 12- and 15-Gyr models). To verify the 
effect of HB stars on Mg b, we have calculated the absorption strength for the following two cases. The 
fitting function of S07 instead of W94 yields, for 12-Gyr SSP with [Fe/H] ~ -1.0, the Mgft sti-ength that is 
decreased by 0.28 A due to the effect of HB stars. The line strengths produced by high-resolution spectra of 



Munari et al.l (|2005r) shows that Mgb decreases by 0.36 A under the same condition. The results imply that 
the effect of HB stars on Mgb, albeit being relatively small compared to H/3, is not negligible and the use 
of Mg b as a direct tracer of metallicity should be with more caution and require modification. It is clear 
that, even after the enhanced a-elements are applied to SSPs, the HB effect on the absorption indices still 
dominates the metal-poor regime irrespective of the [a/Fe] values. 



3.3. Comparison with Other Models 

In this Section, we further demonstrate the characteristics of the YEPS model by comparing it with 
other EPS models (the W94 and TMB03 models). In order to compare the three different models simul- 
taneously, we use the same metallicity scale [Z/H], which is defined as [Z/H] = log(Z/Z0)-log(H/H0), 
where Zq and H© are respectively the total metallicity and Hydrogen content of Sun. Figure |7] displays 
YEPS (solid lines), TMB03 (dashed lines), and W94 (dotted lines) models. Although the three models 
employ heterogeneous sets of stellar evolutionary tracks and EPS modeling schemes, they agree well with 
one another, except for indices CNi, CN2, G4300, H/3, H7, and ¥16. As discussed in previous Sections, the 
YEPS model shows different features in these indices because the indices are particularly sensitive to blue 
HB stars. 

Comparison of the YEPS and TMB03 models for the a-element-enhanced cases shows that the indices 
insensitive to blue HBs are generally in good agreement with each other. In particular, Fe4531, Fe5015, 
(Fe), Fe5406, Fe5709, and Fe5782 agree well. However, the two models do not match for CNi, CN2, 
Ca4227, G4300, Fe4383, Ca4455, and C2466^ The main reasons for the difference are as follows: First, 



^Comparison of the YEPS model with the most recent model of Lee et al. 1 2009^ shows good matches for these indices. When 
the age, [a/Fe], and Z of model are assumed to be 12 Gyr, 0.3, and 0.018, respectively, compared to the scaled solar model, the 
index changes (A/) of CNi, CN2, Ca4227, G4300, Fe4383, Ca4455, and C2 4668 are, respectively,are -0.059, - 0.061, 1.01 -0.52 , 
-2.07, -0.28 and -1.84 for the YEPS model, and -0.054, -0.059, 0.43, -0.78, -2.37, -0.35, and -2.57 for the model of lLee et alj J2009h . 
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the two models use the different prescription of a-elements. The TMB03 model uses fixed C, N, and Ca 
with enhanced a-elements (O, Mg, Na, Ne, S, Si, and Ti) and depressed iron-peak elements (Cr, Mn, Fe, 
Co, Ni, Cu, and Zn). Th e YEPS model, on the other hand, adopts the enhanced ratios of a-elements from 



VandenBerg et all (l2000h . i.e., fixed C and N, enhanced O, Ne, Na, Mg, Si, P, S, CI, Ar, Ca, and Ti, and 



depressed Ar and Mn . The strong Ca4227 of the YEPS, for instance, can be explained by the way we treat 
Ca as a mer nber of the enhanced a-elements group. Second, TMB03 is based on the scaled-solar abundance 
isochrones (|Cassisi et al.lll997l : lBono et al.lll997l : ISalasnich et allbOOOh . The model incorporates the effect 
of a-elements enhancement by depressing iron-peak elements (e.g., Fe and Cr) in a way that satisfies the 
scaling relation [Z/H] = [Fe/H] +0.94[a/Fe]. The effect of the depressed iron-peak elements is not enough 
to mimic the effect of stellar evolution with a depressed iron-peak elements. The YEPS model, on the other 
hand, uses stellar libraries with enhanced a-elements, naturally reproducing the effect of depressed iron- 
peak elements on iron absorption indices without adjusting the iron-peak elements in the stellar atmosphere. 



3.4. Comparison with Observations: Globular Clusters in the Milky Way and M31 



Figures [8] to [TT] present the comparison of the YEPS model with the observed data on the GCs in the 
MW and M31. The model grids in Figures [8] to [12] are for various ages (solid lines) an d metallicities (dotted 
lines). We choose [a/Fel = 0.150 to cons ider the [a/Fe] distribution for the MWGCs jMaraston et allboosl : 
Mendel et al.ll2007l : lwoodlev et al.ll201ol) and the observed mean [a/Fe] of M31 GCs ( ra/Fe1 = 0.14 di 0.04 ; 



Puziaetal 



2005h . Blue triangles represent young (< 1 Gyr) GC candidates in M31 (|Beasley_etalJ[2004|). 



Figures [T3]to[T7]exhibit the line strengths as a function of [Fe/H] for given ages but for differ ent a-elements 
abund a nce ([ a/Fe ] = 0.0, 0.3 an d 0.6). Data for the MW and M31 GCs are obtained from Ischiavon et all 
(|2005l . I2OI2I ') and iBeasley et all ([20041 '). respectively. For t he fair comparison with old age (> 12 Gyr) 
models, we have taken out from young M31 GC candidates ([Beasley et al.ll2004l) in these Figures. With a 
few exceptions, the old GCs in MW and M31 populate well along the 12-Gyr model for the most indices. 
The theoretical predictions of the YEPS model for SSPs at 12 Gyr are given in Tables [5]-(7J 

Notes on individual indices shown in Figures [8] to [T7] are as follows. 

• H/?, H7A, , 116a, and Y16f — Overall shapes of model Balmer lines are in good agreement with 
the observation in Figures [H [3 [131 and[T4l A little offset between metal-rich GCs and the model predictions 
for high order Balmer indices (Figure [8] and [9)1 indicates the [a/Fe] bias between metal-poor and metal-rich 
GCs in the MW and M31 (see Fi gure [T3] and [T4|) . In the metal-poor regime, the 12-Gyr model GCs contain 
hot (> 8,000 K) HB stars that lead to stronger Balmer lines than the model without HB stars. The extensive 



Recent spectr oscopy of individual stars in the MWGCs confi rms that the second generation population in most clusters are 
depleted in oxygen jGratton et al ]|2004ICarretta et alJIlOOsillOOgh . We can expect a similar chemical inhomogeneity in M3 1 GCs 
if the MWGCs are the local counterparts of extragalactic GCs. However, in our models for enhanced a-elements, the abundant 
element oxygen is included in the enhanced a-element group. This would explain why GCs in the MW show better fit with models 
having apparently lower [q/Fel valu es compared to the [a/Fe] measurement based on the [Ca/Fe], [Si/Fe], and [Ti/Fe] of red giant 
stars in the cluster jPritzl et alJuOOST) . 
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study of GCs in M31 by iCaldwell et all (1201 ID confirms that Balmer indices of GCs with blue HBs (see 
their Figure 10) are on average stronger by AH/3 ~ 0.6 A and ARSf ~ 2.0 A. As a consequence, the model 
line is highly inflected, reproducing the observed behavior of the index-index relations. 

The effect of a-elements on Balmer indices differs between H/3 and higher-order Balmer indices. Fig- 
ure [13] shows that H/3 is hardly affected by any o-elements. In contrast, ¥16 and H7 are sensitive to C, Mg, 
Fe, and total metallicity (Figures [T3] and [141). The H5a index is the most sensitive to a-elements enhance- 
ment in the metal-rich regime because of its sensitivity to Fe abundance of giant stars (K05). Note that the 
effect of HB stars on the Balmer indices is universal regardless of a-elements enhancement. 

• Mgi, Mg2, and Mgb — The YEPS model well reproduces the observations of Mgi, Mg2, and Mgb. 
The hnes are simultaneous affected by a-elements, [a/Fe] does not change the overall shape of, for instance, 
the MgZj-Mg2 relation (Figure [141). 

• Fe5270, Fe5335, and (Fe) — The models and observations show a remarkable agreement in Figure[9] 
The lines are widely used as indicators of iron abundance, yet are fairly sensitivity to [a/Fe] (Figure [141). 

• Fe4383 — The YEPS model well reproduces Fe4383 for GCs both in the MW and M31 (Figure 
[9]). Since Fe4383 traces Fe and the total metallicity simultaneously, the Fe4383 model is very sensitive to 
variations of [a/Fe] (Figure [141). 

• Fe4531 — Our model of Fe4531 agrees well with M31 GCs (Figure [TOl). This index is a good iron 
abundance indicator that is relatively less sensitive to [a/Fe] (Figure [T5l). 

• Fe5015 — The YEPS model shows a good agreement with the M31 GCs (Figure [TOt. The Fe5015 
strength is generally involved with the Ti and Mg abundance in low MS stars and giant stars, as well as total 
metallicity near TOs. Yet, the effect of [a/Fe] on Fe5015 is only modest (Figure [T5l). 

• Fe5406, Fe5709, and Fe5782 — GCs in the MW and M31 show good matches with the models. The 
slight offset of MWGCs in Fe5782 (Figure [TT]) is likely due to the well-known problem of miscalibration 
of fitting functions (TMB03). The a-element sensitivity of these indices are similar to those of Fe5270 and 
Fe5335 (Figures[l5]and[l6l). 

• CNi and CN2 — The models for CNi and CN2 do not show good fits to the observational data in 
Figure [TT] Compared to observational data, our models predict on average 0. 1 mag lower values in the 
metal-rich regime. The in ferred [a/Fe] seems u nreasonably high and likely due to the poor calibration of 



fitting functions (TMB03: lLee & Wortheyll2005h . The poor calibration of fitting functions can be explained 



by the effect of CN ONa anticorrelation found in the MWGCs with the second generation stellar population 



([Coelho et al.[[201 ih . CNi and CN2 indices, even for the a-element enhanced models, are affected fairly by 



blue HB stars (see sohd and dashed lines in Figure [T6l). 

• Ca4227 — Our models for Ca4227 show a small offset with observations. Observatio ns occupy 



slightly low EWs compared to the model in Figure[TT]and Figure[T6] As suggested by TMB03 and lLee et al. 
the enhanced C and N model can depress EW of Ca4227 but we do not consider the case of the en- 



hanced C and N model in this study. Note that more realistic model with the effect of CNONa anticorrelation 
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would have decreased strengths of Ca4227 dCoelho et alj|201 ih and this may explain the offset between the 



YEPS model and observations. Since Ca4227 is very sensitive to Ca and C, but is not so to Mg, the Ca4227 
line deepens with increasing [a/Fe] (Figure [TQ. 

• Ca4455 — The Ca4455 model for [a/Fe] = 0.15 predicts on average 0.3 A higher EWs compared 
to the observations (Figure [T2l). and the observations are better reproduced by the model with [a/Fe] ~ 0.3 
(Figure [TtI). Since Ca4455 is insensitive to the variation of any of elements studied in K05, the effect of 
[a/Fe] shown in Figure [17] come solely from the variation of Fe due to the increased [a/Fe] at a given Z. 
This explains why Ca4455 shows a different response compared to Ca4227 even though we treat Ca as a 
member of enhanced a-elements group (Table |2l). 

• G4300 — The YEPS model with enhanced a-elements ([a/Fe] = 0.15) for the G4300 shows a httle 
offset from the observations (Figure [12)) . Rather, GCs in the MW and M31 follow the model with [a/Fe] 
= 0.0 well (Figure [TT]). Given that the other lines agree with the observations, our G4300 model seems too 
sensitive to [a/Fe]. 

• C24668 — The YEPS model for the C24668 index offsets from the observations (Figure [12]) . The 
model with enhanced [a/Fe] (Figure [TT]) shows great sensitivity to [a/Fe]. C24668 is a measurement of C 
abundance, and is insensitive to all other a-elements. Because the YEPS model with enhanced [a/Fe] used 
fixed C abundance, the strong sensitivity to the a-elements enhancement in Figure [TT] comes mainly from 
the [Fe/H] variation due to the a-element contents at a given Z. 

• NaD — The YEPS model for NaD index is weaker than the observations. Part of the reason for this 



discrepancy is Na absorption by the interstellar medium (TMB03). As suggested by iCoeUio et all (120111) . 
more realistic model with increased Na abundance caused by the observed CNONa anticorrelation in the 
second generation population would increase the strength of NaD. The Na sensitivity of NaD plays an 
important role in the a-element enhanced model by decreasing the strengths of NaD at given Mg b (Figure 

[n]). 

Based on the comparison shown in Figures [8] to [TT] we identify several Lick indices that are most 
appropriate for the estimation of metallicity, age, and a-elements enhancement of stellar populations. (1) 
For determination of iron abundance, Fe4383, Fe4531, Fe5015, Fe5270, and Fe5335 are recommended; (2) 
For age-dating, Balmer indices such as H/3, H7A, H777, HJ^, and HJ^ ; and (3) For measuring a-elements 
enhancement, Mgi, Mg2, and Mg b. 



4. DISCUSSION 

4.1. Estimation of Age and Metallicity with YEPS 

In this Section, we discuss how the effects of HB stars and a-elements enhancement in the model are 
combined to affect the age and metallicity estimation of SSPs. To this aim, we select three typical MWGC , 
47 Tucanae (NGC 104), NGC 6284, and M67 (NGC 2682), representing old metal-rich, old metal-poor, and 
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young metal-rich GC populations, respectively. These clusters have well-studied CMDs and show no strong 



evidence of multiple stellar populations reported by recent studies (ILee et al 



I997I : Lee et alJl2009l : ban et alJl2009l : Iperraro et alJboogI : Ipiotto et al 



2007 



1999; 



Lavden & Saraiedini 



Moretti et al 



200' 



i 



Figure [18] displays the observed and synthetic CMDs for 47 Tuc, NGC 6284, and M67. The free pa- 
rameters used to match the synthetic CMDs with the observed ones are age and metallicity. To simulate 
observational errors in our model, we also carried out Monte Carlo simulations based on the actual observa- 
tional uncertainties. The best-fit parameters of each model are summarized in Table |4l 

The left column of Figure [19] compares the absorption indices between the observations and models for 
the typical GCs in the MW. Since the Y EPS models wit hout HB stars generate very similar results to other 
models with red clump single-mass HB dLee et al.ll200d '). we present the YEPS model with and without HB 
stars in Figure [19] (1) to highlight the effect of HB, and (2) to compare our model with other models with 
red clump single-mass HBs. The top panel shows that the H/3 indices of the GCs are reproduced better 
by the model with the HB effect and the ages based on the absorption lines are in better agreement with 
the ages (Table |4]) derived from synthetic CMDs in Figure [18] In particular, the systematic variation of 
HB morphologies with respect to metallicity and age is essential to explain the enhanced H/3 strength of 
NGC 6284 — an old, metal-poor system with well-developed blue HB stars. The age of NGC 6284 would be 
estimated to be 8-9 Gyr without the HB effect, which is inconsistent with its derived age (13.1 Gyr) from 
the MSTO and HB morphology in Figure[T8] On the other hand, the absorption strengths of 47 Tuc and M67 
show reasonable agreement with the YEPS models for 12-Gyr (solid cyan lines) and 3.5-Gyr (between solid 
red and orange lines) GCs, respectively. The models for both metal-rich and young clusters possess red HB 
stars. For a 3-Gyr model GC, the red HB reduces H/3 by 0.2 A, which corresponds to ~ 1 Gyr. Therefore, in 
order to derive accurate ages of stellar populations based on the EPS model prediction of B aimer strengths, 
one should check first whether the EPS model has well calibrated HBs for both blue-HBs (e.g., NGC 6284) 
and red-HBs (e.g., 47 Tuc and M67). 

The righ t column of Figure [T9] is similar to the top panel of the left column, but shows the age dating 
of M31 GCs ( Beaslev et al.ll2004h . In these panels, based on the a-elements enhancement of GCs in M31 
(|Puzia et al.ll2005h . we have used the model of [a/Fe] = 0.15 for the comparison. The metal-poor GCs in 
M31 exhibit stronger H/3 by ~1 A than the metal-rich counterparts. The upper panel of the right column 
shows, when the effect of HBs is not considered, the models give ~ 5 Gyr younger ages for the metal-poor 
GCs compared to the metal-rich GCs in M31. With the general trend of age-metallicity relation in mind, it 
would be difficult to interpret that the mean age of metal-poor GCs is younger than that of metal-rich GCs. 
By contrast, the lower panel shows that the systematically enhanced Balmer indices of the metal-poor M3 1 
GCs are well reproduced by our single age, 12-Gyr model GCs with HB stars. Indeed, Perina et al. (2009) 
directly detected well developed blue HB stars in the metal-poor GCs in M31 (e.g., BO 10, B220, B224, and 
B366) based on the HST/ACS CMDs. 

Back in the left column of Figure [19] the middle and bottom panels highlight the effects of HBs and 
a-elements enhancement on the metallicity determination of SSPs. The effects of a-elements and HB stars 
go in the opposite direction on absorption strengths of Mgb for given [Fe/H]'s; the inclusion of blue HB 
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stars decreases indices, whereas the enhancement of a-elements increases indices. While the HBs exert only 
a marginal effect on determining [Fe/H], the a-elements have a marked effect. For instance, the effect of 
enhanced a-elements ([a/Fe] = 0.3) on Mgb is about 10 times greater than that of HBs when the age and 
[Fe/H] are assumed to be 12 Gyr and 0.0 dex, respectively. Note that the model line of Mg bfor 12 Gyr GCs 
with [a/Fe] = 0.3 is almost identical to that for 4 Gyr GCs with [a/Fe] = 0.0 due to strong sensitivity of the 
Mgb to the a-elements. In order to determine accurate [a/Fe], therefore, the age dating and the metallicity 
determination of SSPs should be carried out at the same time. The (Fe) index, on the other hand, is less 
sensitive to [a/Fe] than Mg b, which makes (Fe) a more accurate [Fe/H] indicator for the SSPs than other 
indices sensitive to [a/Fe] (e.g., Mgb and Mg2). 

Since the age dating with Balmer indices is usually applied to the samples of elliptical galaxies (e.g.. 



aince tne age aatmg witn maimer indices is usually applied to tne samples or elliptical galaxies (e.g., 
Trager et allboool i lcaldwell et al.l2003l : lThomas et alboOsl i lNelan et albooilrrager et allboOsl i lGraves etal 
2009h . the effect of HB stars is also important for composite stellar populations. The effect of HB stars on 



the integrated observables of giant elliptical galaxies should be limited because their mean metallicities are 
as high as [Z/H] > 0.0, for which hot HB stars are rare. But the metallicity spread of ellipticals, nevertheless, 
allows them to have a cert ain fraction of low-metallicity stars and thus hot HB stars accordingly ([Park & Lee 
1997l: lchunget"^ l201 1|). Moreover, the effect of hot HB stars in dwarf elliptical galaxies is more important 
because the metallicity o f dwarf ellipticals a re generally lower than that of giant ellipticals. In the context of 
galaxy downsizing (e.g.. lCowie et al.lll99q) . the age dating of giant and dwarf elliptical galaxies with well 
calibrated HB models is crucial issue for determining the formation history of elliptical galaxies. We will 
fully discuss this issue in our forthcoming paper (C. Chung et al. 2012, in preparation). 



4.2. Distributions of Absorption Indices in Extragalactic GCs 

Ever since the recogn i tion of bimodal broadband color distributions of GCs in massive early-type galax- 



ies (2 


:eDf & Ashmanlll993 


:Geisleret all 1996: 


Forbes et al. 


1997: 


Gebhardllll999: 


KundJboOl: 


Larsen et al. 


2001 


; Jordan et al. 


2002:1 


'uzia etal.l2004: 


West et al.l2004l: 


Strader et albood: 


Pens et al.l2006 


: Harris et al. 


2006 


; Jordan et al. 


2OO9I: 


Liu et al. 


2011). the phenomenon has been interpreted as the presence of two GC 



subsystems within i ndividual galaxies. T hree major ideas have been put forw ard to explain it, in cluding the 



merg er model (e.g.. Ashman et al. 1992 ). the in-situ (e.g., Forbes et al. 19981) and accretion (e.g.. Cote et al. 
199 w scenarios. 

More recently, however, Yoon et al. (2006) and Yoon et al. ( 2011al lbl) suggest an alternative explanation 
for the GC color bimodality that does not necessarily invoke two GC subpopulations. Yoon et al. show that 
the theoretical metallicity-color relations are inflected and that such relati ons can generate bimodal color 
distributions from broad metallicity spreads, even if they are unimodal. The I Yoon et al.l (120061) model of GC 
colors indicates that the HB effect is the most important for the inflection on the metallicity-color relations. 
HBs also have strong effect on the absorption index versus metallicity relations (IMRs), and thus one can 
put the Yoon et al. explanation to the test by examining the metalhcity-index nonhnearity and the resulting 
index distributions. 
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4.2.1. Conversions from Metallicity Spreads to Index Distributions 



Figure [20] shows the conversion from metallicity spreads to index distributions via the IMRs. Three 
different ages are selected that show little (5 Gyr) and significant (12 and 13 Gyr) inflection on the IMRs. 
We assume the underlying metallicity distribution functions (MDFs) of 10^ model GCs to be box-shaped 
and perform the Monte Carlo simulations for index distributions. The simple MDFs should allow us to see 
the pure effect of the IMR projection. The indices of each model GC are calcula ted from its [Fe/HI v alue 
via the corresponding IMRs and then typical errors estimated from observation of iBeasley et al.l (|2004h are 
randomly added. Our simulations with the non-inflected IMRs do not produce bimodal index distributions. 
For example, Fe4531 and Fe5782 of the 5 Gyr models are trapezoid-shaped. In many cases, the index 
histograms have a sharp peak with a long tail as the IMRs are broken roughly into parts — the metal-poor, 
steep section and the metal-rich, shallow section. On the other hand, the index distributions produced by 
the highly inflected IMRs clearly show bimodality. For example, the IMRs for G4300, H/3, H7a,f, and 
H(5A,f of the old GCs have a very shallow slope at [Fe/H] ~ -1.0 between two steeper slopes. The inflection 
brings about a dip on the index distributions by projecting equidistant metallicity intervals onto larger index 
intervals. 

Figure EUrepeats a similar experiment but uses Gaussian MDFs. We test for the two Gaussian MDFs of 
([Fe/H]) = -1 .0 and -0.7 with the same dispersion of (T[fg///] = 0.55. In this simulation, we also have applied 
10^ model GCs for the given Gaussian MDFs. Even with Gaussian MDFs, the projected index distributions 
are double-peaked for G4300, H/3, H7A,f , and of 12-Gvr mode l GCs. When a MDF with ([Fe/H]) 

= -1.0 is used, for example, the KMM algorithm (lAshman et al.lll994)) strongly prefer two peaks for these 
histograms with p- value ~ 0.0. For the projected H/3 distribution under the assumption of 12 Gyr model, 
the two peaks are located at 1.681 A and 2.535 A with a number fraction of 48.6 % and 51.4 %, respectively. 

It has been claimed that the bimodal distributions of metal-line indices (e.g., MgZ? and [MgFe]') of 
GCs in ea rly-type galax i es are the evidence that GCs have two subsystems with different metallicities 
(ICohen et al. 1998. 12OO3I ; Istrader et allboO?! ; Iwoodley et allboiol) . Interestingly, however, even with uni- 
modal MDF, the index distribution of Mg b, frequently used as a metallicity indicator, also shows a bimodal 
index distribution that is consistent with the observations. Our KMM test for the projected Mg b distribution 
of the 12- and 13-Gyr models supports bimodal distributions with the p-value of 0.0. The test implies that 
the HBs exert an appreciable effect not only on the Balmer lines but also on the metallicity-sensitive lines. 
Without assuming two subpopulations in GCs, the projection effect can reproduce the observed bimodal 
Mg b histograms. 



4.2.2. Absorption Index Distributions ofMBl Globular Clusters 



In order to compare the observational index distributions with the simulated ones, we choose M31, 
the nearest large galaxy, which has spectroscopic su rveys of a great number of GCs with reasonably small 
observation errors. M31 GCs ([Beasley et al.l 120041) are of 18 and 33 % smaller errors respectively in H/3 
and (Fe) at given absolute magnitudes compared to the M87 GC spectroscopy with Subaru (S. Kim et al. 
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2012 in preparation). We excluded young cluster candidates in lBeasley et alj (120041) to avoid young cluster 
contaminations in the MDF (see blue triangles in Figure [8]to[T2l). 

Figure [22] shows the index-index diagrams for (Fe) and Balmer indices. The observed index distribu- 
tions of M31 GCs and the projection simulations based on the YEPS and TMB03 IMRs are shown as the 
histograms along the y- and A;-axes, respectively. We, again, make a simple assumption of Gaussian MDFs 
with ([Fe/H]) = -1.0 and o"[pe/H] - 0.55. As explained in detail in Section 4.2.1, the principle of the MDF 
projection with IMRs is to transfer SSPs from the metallicity space to the index space. Red circles with an 
uniform metallicity intervals (A[Fe/H] ~ 0.3) on the YEPS IMRs in the bottom panel show this principle of 
projection. All Balmer indices show large index space between [Fe/H] = -1.1 to -0.5 and the corresponding 
index distributions also show relatively low number frequencies. However, the (Fe) index, which keeps rel- 
atively constant index intervals at fixed [Fe/H] intervals, does not show a significant change in the projected 
distribution. As a result, the projection simulation with the YEPS model reproduce the unimodal (Fe) and 
bimodal Balmer line distributions at the same time using a single Gaussian MDF. We stress that if the obser- 
vational errors are small enough then the (Fe) distributions can also show a weak bimodal feature. Indeed, 
the rece nt observation of GCs in M31 by lCaldwell et all (|201l[ ) show a weak bimodality in the distribution 
of (Fe) dKim et al.|[2012l) . The TMB03 model, on the other hand, gives almost straight IMRs and does not 
reproduce the observed distributions of the (Fe) and Balmer indices, simultaneously. No matter how small 
the errors are, the single Gaussian MDF projection based on the models without systematic HB effect can 
not reproduce the observed distributions for both (Fe) and Balmer indices. 

Figure [23] shows the index-index diagrams for Mg h and Balmer indices. The observed and simulated 
Mgft and Balmer distributions are displayed along the y- and x-axes, respectively. As discussed in Section 
4.2.1, Mg h is fairly affected by blue HBs and thus the distribution of MgZ? shows weak bimodality. Many 
other observations of GCs in early type galaxies also show simi lar results and the bimodal Mg h distributions 
are often interpreted as the evidence of the bimodal GC MDFs Jcohen et alii 9981 l20oi : lstrader et allboO?! : 
Woodley et al.[[2010l) . The comparison of the YEPS model to M31 GCs demonstrates that the HB effect on 
Mg h is certainly non-negligible. Greater caution, therefore, is required in deriving GC metallicity directly 
from Mg h. 

The spectroscopic line indices, as more detailed probes of stellar populations than broadband colors, 
contain abundant information on the structure of GC systems. Our simulations targeting at the old GCs in 
M31 suggest that two distinct groups found in Balmer and Mg h indices can be due simply to the inflected, 
nonlinear IMRs. The result favors a unimodal [Fe/ H] distribution of M31 old GCs, in line with the GC 



Yoon et al 



J2006l . l2011al lbl). Absorption indices 



are 



structure of extragalactic GC systems suggested by 
generally subject to larger observational uncertainties and spectroscopic samples are still small compared 
to photometric samples with broadband colors. More precise spectroscopic observations of greater number 
of extragalactic GCs by next-generation telescope projects, such as the Giant Magellan Telescope (GMT), 
Thirty Meter Telescope (TMT) and European Extremely Large Telescope (E-ELT), are highly anticipated in 
this regard. 
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5. SUMMARY 

We have presented an updated and blown-up version of the Yonsei Evolutionary Populations Synthesis 
(YEPS) model for spectroscopic absorption indices of simple stellar populations. The characteristics of the 
YEPS and its appUcations are summarized as follows. 

1 . The YEPS has included detailed HB models, which reproduce the observed HB morphology of the 
Milky Way GCs and varies with respect to metalUcity, age, and a-elements enhancement. 

2. The YEPS has incorporated the a-element variation by using the Y^ stellar library with enhanced 
a-elements and the response functions of a-elements by K05. 

3. The YEPS is in good agreement with other EPS models, except for the indices sensitive to hot stars 
(CNi, CN2, G4300, H/3, H7, and H5). 

4. The YEPS reproduces well the observed absorption features of GCs in the MW and M31, including 
the strengthened Bahner absorptions by the effect of blue HBs. 

5. When the observed variation of HB morphology with metallicity in the MWGCs is included in 
the models, Balmer indices of SSPs do not monotonically decrease with increasing metallicity at given age 
because of blue HB stars in the metal-poor regime. Therefore, the age dating of old stellar systems based on 
Balmer indices suffers from age degeneracy in the metal-poor regime. 

6. The contribution of HBs to absorption strengths is not limited to Balmer indices but influences the 
indices pertinent to iron and a-elements. As a consequence, the most index-metalUcity relations of YEPS 
are inflected and nonlinear. 

7. We have simulated, for the first time, the index distributions of GCs using the YEPS index- 
metallicity relations. The distributions of Balmer and Mgb indices constructed under the unimodal MDF 
assumption show clear bimodality, which can be viewed as a close analogy to the well-known bimodality 
found in the broadband color distributions of extragalactic GC systems. 
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Fig. 1. — The variation of HB morphology as functions of metallicity and age . Filled and open circles and 
triangles are GCs in three radial zones of the MiUcy Way ( Lee et all 1994 . 2007 ). Solid blue lines are theoret- 
ical isochrones from 10, 11, and 12 Gyr YEPS model. Blue and red arrows indicate the selected metallicities 
([Fe/H] = -2.5, -1.7, -0.9, and -0.3) for synthetic color-magnitude diagrams (CMDs) in Figured 



-25- 



n — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

- 10 Gyr, [Fe/H] = -0.3 
HBtype = - 1.000 




I I I I I I I I I I I I I n I I 

- 10 Gyr, [Fe/H] = -0.9 
HBtype = - 1.000 




I I I I I I I I I I I I I I I 
-10 Gyr, [Fe/H] = -1.7 
HB type = - 0.688 




I I I I I I I I I I I M I I I I 
-10 Gyr, [Fe/H] = -2.5 
HB type = 0.833 




J I I I l_ 



"1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

-12 Gyr, [Fe/H] = -0.3 
HBtype = - 1.000 




I I I I I I I I I I I I I n I I 

- 12 Gyr, [Fe/H] = -0.9 
HB type = - 0.982 




I I I I I I I I I I I I n i I I I 

-12 Gyr, [Fe/H] = -1.7 
HB type = 0.932 




I I I I I I I I I I I M I I I I 
-12 Gyr, [Fe/H] = -2.5 
HB type = 0.989 




J I I I l_ 



-0.5 



0.5 
B-V 



-0.5 



0.5 
B-V 



Fig. 2. — The effect of metallicity and age on the morphology of HB stars in the CMD. IVIetallicity of 
each CMD increases from the bottom to the top panels. CMDs on the left and right sides are 10 and 12 
Gyr models, respectively. Solid lines are isochrones and dots are corresponding HB stars at given ages and 
metallicities. 
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Fig. 3. — H/3 fitting functions witli respect to the gravity and temperature (0 = 5040/Teff). Black dots are 
H/3 measurements of 460 stellar spectra used for the construction of W94's fitting functions. Left panels are 
fitting functions for giant stars (with logg < 3.6) and right panels are fitting functions for dwarf stars (with 
logg > 3.6). Upper and bottom panels are fitting functions of W94 and S07, respectively. Red, green, and 
blue solid lines indicate metallicities of [Fe/H] = -1.0, -0.3, and 0.1, respectively. 
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Fig. 4. — The effect of HB stars on the absorption indices of CNi, G4300, H'Jf, and H/3. Red, green, and 
blue lines indicate, respectively, SSP models with ages of 10, 12, and 14 Gyr. Solid and dashed lines are 
YEPS models with and without HB stars, respectively. The models are for [a/Fe] = 0.3. 
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Fig. 5. — The effect of various input parameters (stellar libraries and fitting functions) on the strength of 
H/3. Three color codes indicate different ages of 10, 12, and 14 Gyr. Solid and dashed lines are YEPS 
models with and without HB stars. Models in upper panels are calculated based on fitting functions of S07 
and models in bottom panels are con structed using fit ting functions of W94. Stellar evoluti onary tracks used 
in the left and right panels are (IKim et al.l l2002h and BaSTI (IPietrinfemi et all 120041) stellar libraries, 
respectively. 
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Fig. 6. — The effect of a-elements on the HB type, H/3, and Mg b indices. Green, red, and blue colors are 
the YEPS models for [a/Fe] = 0.0, 0.3, and 0.6, respectively. The left panel is the same plot as Figure[T]but 
for various a-element enhancements. Solid, dotted, and dot dashed lines in the left panel are HB types for 3 
different ages of 15, 12, and 9 Gyr, respectively. Middle panels are the H/3 models with and without HB stars 
for different ages. Solid and dashed lines in the middle and right panels represent the YEPS models with and 
without HB stars. From top to bottom panels in the middle and right panels, the ages of the YEPS models 
are 9, 12, and 15 Gyr. Right panels are the YEPS models for the Vlgb index with the same parameters as in 
the middle panels. 
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Fig. 7. — Comparison of the YEPS models with W94 and TMB03 models. All models are the same age of 
12 Gyr. Solid lines with green, red, and blue colors are predictions of the YEPS model for [a/Fe] = 0.0, 0.3, 
and 0.6, respectively. Green, red, and blue dashed lines represent TMB03 models for [a/Fe] = 0.0, 0.3, and 
0.5. Dotted lines in green indicate W94 models with scaled solar abundance. 
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Fig. 8.— Comparisons of the YEPS model ([a/Fe] = 0.15) with GCs in the MW and M31. Absorption 
indices of H/3, H'Ja, ^"/f, and ¥16a are presented in this figure. Solid lines are the same age SSP model 
and dotted black lines indicate iso-metallicities from [Fe/H] = -2.5 to 0.5. Each solid line corresponds to 
the age between 1 to 15 Gyr in steps of 1 Gyr from red to blue colors. Thick lines with turquoise and dark 
blue colors are t he YEPS model for the age of 12 and 13 Gyr, respectively. Red squares a re GCs in the 
MW o bserved by lSchiavon et al.l (|2005l . l2012h . Black and blue triangle s are GCs in M31 fro m lBeasley et al. 
(|2004h . Blue triangles are young cluster candidates of M31 based on iBeasley et al.l (|2004l '). A little offset 
between metal rich GCs and the model predictions for high order Balmer indices indicates the [a/Fe] bias 
between metal-poor and metal-rich GCs in the MW and M31. 
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Fig. 13. — Comparisons between the YEPS model with enhanced a-elements and GCs in the MW and 
M31. Absorption indices of H/3, H7A, Hjf, and Y16a are displayed in this figure. Green, red, and blue lines 
represent the YEPS model for [a/Fe] = 0.0, 0.3 and 0.6, respectively. Solid and dashed hnes are the YEPS 
models for 12 and 14 Gyr, respectively. Observed GCs are from the same data used in Figures [8] to [T2l but 
we have excluded young M31 GC candidates (blue triangles in Figures from[8]to[T2l) for the fair comparison 
with old age (> 12 Gyr) models. 
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Fig. 18.— Comparison of the observed 47 Tuc, NGC 6284 (IPiotto et all 120021) . and M67 (IStassun et al. 
2002h with the synthetic CMD models (bottom panels). Observed CMDs are displayed with matched Y^- 
isochrones. 
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Fig. 19. — Comparison of the YEPS model with the observed absorption indi ces of NGC 6284, 47 Tuc 



dSchiavon et al.ll2005l . l2012r) . and M67 (ISchiavonll2007l) . as well as GCs in M31 (IBeaslev et al.ll2004h . Blue, 
red, and cyan squares in the left panels are observed as integrated absorption indices of NGC 6284, 47 Tuc, 
and M67. Solid and long dashed lines in all panels are models with and without HB stars. Red, orange, grey, 
green, turquoise, and blue lines in the left top panel represent the H/3 models for the ages of 3, 4, 5, 8, 12 and 
13 Gyr with [a/Fe] = 0.3. Grey and black lines in the middle and bottom panels on the left side are models 
for [a/Fe] = 0.0 and 0.3 at the age of 12 Gyr. Light blue and blue lines in the same panels are models for 
[a/Fe]=0.0 and 0.3 at the age of 3 Gyr. In the right p anels, the effect of HB stars on the age-dating with the 
H/3 is displayed with M31 GCs (IBeaslev et al.ll2004t) . Model lines in the right panels are for [a/Fe]=0.l5. 
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Fig. 20. — The test of the projection effect with the box-shaped metalUcity distribution function (MDF) on various 
Lick absorption indices. The box-shaped MDFs of 10'' model GCs are shown along with y-axis. Black dashed, solid 
and dotted lines in each plot indicate 5, 12, and 13 Gyr models with [a/Fe] = 0. 15 for given Lick indices. Red dashed, 
solid, and dotted lines in the bottom of each plot show the result of box-MDF projection to various Lick absorption 
indices. Absorption index distributions of each line type are corresponding to the projection through the models with 
the same line types. 
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Fig. 21. — Same as Figure |20] but for the projection tests with the single Gaussian MDFs of ([Fe/H]) = 
-0.7 (red) and -1.0 (blue) on various Lick absorption indices. The 10^ model GCs with single Gaussian 
MDFs are projected on various absorption indices via IMRs of the YEPS model. 
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Fig. 22. — Comparison between index distributions of M31 GCs and simulated index distributions. Upper 
panels are observed histograms of M31 GCs in H/3, ll'yA,F, H^a.f , and (Fe) and the YEPS models (red lines) 
for each index are overlaid with GCs in M31. Lower panels are simulated index distributions with the YEPS 
and TMB03 models. Black sohd and dashed Unes represent the YEPS and TMB03 model for [a/Fe]=0.15, 
respectively. Blue solid and dashed lines are simulated index distributions based on the YEPS and TMB03 
model, respectively. Red circles on the solid lines indicate metallicities from [Fe/H] = -2.6 to -0.5 with 
a fixed interval A[Fe/H] = 0.3. Three values of [Fe/H] = -1.0, -0.8, and -0.5 are indicated in the YEPS 
model for H/3. 
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Table 1 . Model Input Parameters 



Input Ingredients & Parameters 



Standard Model 



Comparison Model 



Stellar Library 

Spectral Library 
Empirical Fitting Functions 



Response Functions of a-elements 

Initial Mass Function 

o-elements enhancement, [a/Fe] 

HB mass dispersion, cfm{Mq) 

Reimers (1977)'s mass-loss efficiency parameter, ri 

Assumption of the age of inner-halo MWGCs 



Y stellar libraries 
(Kim et al. 2002 : Lee et al. 2012 in prep.) 
BaSel 3.1 rWestera et al. 20021 
Worthev (1994') 
Worthev & Ottaviani (1997) 
Schiavon (2007) 
Kom et al. (2005) 
Salpeter (x= 1.35) 
0.0, 0.3, 0.6 
0.015 
0.63 
12Gyr 



BaSTI 
(Pietrinfemi et al. 2004) 
BaSel 3.1 (Westera et al. 2002) 
Worthev (1994) 
Worthev & Ottaviani (1997) 
Schiavon (2007) 
Kom et al. (2005) 
Salpeter (x=\3$) 
0.0, 0.4 
0.015 
0.40^ 
12Gyr 
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Table 2: The a-element enhanced patterns of stellar libraries 



a-element 


Grevesse 


y -Isochrones 


BaSTI 




[a/Fe] = 0.0 


[a/Fe] = 0.3 


[a/Fe] = 0.4 


C 


8.55 


8.55 


8.55 


N 


7.97 


7.97 


7.97 


O 


8.87 


9.17 


9.37 


JNe 




0.38 


O TT 

8.37 


Na 


6.33 


6.63 


6.33 


Mg 


7.58 


7.88 


7.98 


Al 


6.47 


6.17 


6.47 


Si 


7.55 


7.85 


7.85 


P 


5.45 


5.75 


5.45 


S 


7.21 


7.51 


7.54 


CI 


5.50 


5.80 


5.50 


At 


6.52 


6.82 


6.52 


K 


5.12 


5.12 


5.12 


Ca 


6.36 


6.66 


6.86 


Ti 


5.02 


5.32 


5.65 


Cr 


5.67 


5.67 


5.68 


Mn 


5.39 


5.24 


5.39 


Fe 


7.50 


7.50 


7.50 


Ni 


6.25 


6.25 


6.29 



Note. — The abundance of elements is listed in logarithmic scale log Nei/Nn + 12. 
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Table 3: Minimum Lick indices based on W94 fitting functions 

Lick Index stellar library BaSTI 



CNi 


-0.329 


-0.301 


CN2 


-0.309 


-0.264 


Ca4227 


-0.654 


-0.642 


G4300 


-5.712 


-5.712 


Fe4383 


^.371 


^.377 


Ca4455 


-0.400 


-0.404 


Fe4531 


-1.448 


-1.417 


Fe4668 


-6.238 


-2.156 


H/3 


-1.726 


-1.726 


Fe5015 


-0.838 


-0.890 


Mgi 


-0.170 


-0.170 


Mg2 


-0.082 


-0.082 


Mgb 


-1.442 


-1.444 


Fe5270 


-2.375 


-2.350 


Fe5335 


-0.275 


-0.294 


Fe5406 


-0.952 


-0.951 


Fe5709 


-1.976 


-1.976 


Fe5782 


-0.805 


-0.804 


NaD 


0.000 


0.000 


TiOi 


-0.068 


-0.018 


Ti02 


-0.055 


-0.008 


H7A 


-12.269 


-12.152 


H7F 


^.292 


^.278 


USa 


-9.384 


-9.139 




-2.359 


-2.249 



Table 4. Input Parameters of CMD model 



Parameters 47 Tuc NGC 6284 M67 



Initial Mass Function Salpeter (a' = 1.35) Salpeter (jc = 1 .35) Salpeter (a = 1.35) 

a-elements enhancement, [a/Fe] 0.3 0.3 0.3 

HB mass dispersion, o-mCMq) 0.015 0.015 0.015 

Reimers ^97T)'s mass-loss efficiency parameter, r] 0.63 0.63 0.63 

Distance modulus, (V-Mv) (mag) 13.35 16.70 9.95 

Galactic reddening, £(B-V) (mag) 0.01 0.32 0.04 

Metal abundance, [Fe/H] -0.73 -1.50 0.00 

Absolute age, f (Gyr) 11.9 13.1 3.5 



Table 5. LICK absorption indices of YEPS simple stellar population model for [Q;/Fe]=0.0 (fitting 

functions of W94 and W97). 



Age = 12.0 



[Fe/H] 


CNi 


CN2 


Ca4227 


G4300 


Fe4383 


Ca4455 


Fe4531 


C24668 


H/J 


Fe5015 


Mgi 


Mg2 


-2.50 


-0.132 


-0.089 


0.213 


0.033 


0.451 


0.166 


0.258 


0.156 


2.672 


0.798 


0.007 


0.021 


-2.40 


-0.126 


-0.082 


0.230 


0.207 


0.426 


0.181 


0.378 


0.010 


2.673 


0.885 


0.007 


0.023 


-2.30 


-0.121 


-0.076 


0.242 


0.342 


0.401 


0.196 


0.490 


-0.100 


2.693 


0.982 


0.006 


0.024 


-2.20 


-0.118 


-0.072 


0.256 


0.465 


0.393 


0.214 


0.600 


-0.202 


2.719 


1.086 


0.006 


0.026 


-2.10 


-0.116 


-0.069 


0.265 


0.563 


0.390 


0.232 


0.703 


-0.285 


2.758 


1.188 


0.006 


0.028 


-2.00 


-0.114 


-0.065 


0.281 


0.688 


0.421 


0.254 


0.813 


-0.359 


2.775 


1.292 


0.007 


0.032 


-1.90 


-0.112 


-0.063 


0.298 


0.810 


0.468 


0.278 


0.922 


-0.417 


2.794 


1.399 


0.008 


0.035 


-1.80 


-0.110 


-0.060 


0.320 


0.958 


0.543 


0.308 


1.035 


-0.454 


2.792 


1.514 


0.010 


0.040 


-1.70 


-0.107 


-0.058 


0.344 


1.104 


0.637 


0.341 


1.149 


-0.465 


2.793 


1.633 


0.012 


0.046 


-1.60 


-0.104 


-0.055 


0.378 


1.301 


0.768 


0.382 


1.273 


-0.453 


2.765 


1.771 


0.015 


0.053 


-1.50 


-0.099 


-0.052 


0.418 


1.538 


0.931 


0.431 


1.405 


-0.406 


2.712 


1.926 


0.018 


0.061 


-1.40 


-0.094 


-0.049 


0.467 


1.813 


1.119 


0.489 


1.545 


-0.312 


2.630 


2.103 


0.021 


0.070 


-1.30 


-0.089 


-0.045 


0.517 


2.080 


1.320 


0.550 


1.683 


-0.179 


2.556 


2.292 


0.026 


0.081 


-1.20 


-0.084 


-0.042 


0.566 


2.339 


1.527 


0.613 


1.816 


-0.018 


2.499 


2.487 


0.030 


0.092 


-1.10 


-0.076 


-0.037 


0.634 


2.718 


1.799 


0.695 


1.976 


0.202 


2.376 


2.719 


0.035 


0.104 


-1.00 


-0.066 


-0.030 


0.705 


3.147 


2.094 


0.785 


2.140 


0.468 


2.226 


2.962 


0.040 


0.117 


-0.90 


-0.055 


-0.022 


0.787 


3.606 


2.448 


0.886 


2.313 


0.799 


2.068 


3.227 


0.046 


0.131 


-0.80 


-0.045 


-0.015 


0.868 


4.028 


2.844 


0.989 


2.479 


1.176 


1.939 


3.492 


0.053 


0.145 


-0.70 


-0.035 


-0.007 


0.959 


4.417 


3.277 


1.099 


2.648 


1.613 


1.815 


3.773 


0.060 


0.158 


-0.60 


-0.028 


-0.001 


1.040 


4.695 


3.658 


1.200 


2.799 


2.054 


1.730 


4.038 


0.068 


0.172 


-0.50 


-0.018 


0.008 


1.108 


4.946 


3.952 


1.295 


2.932 


2.463 


1.660 


4.332 


0.074 


0.185 


-0.40 


-0.003 


0.023 


1.182 


5.171 


4.269 


1.394 


3.069 


2.940 


1.604 


4.634 


0.079 


0.199 


-0.30 


0.009 


0.036 


1.286 


5.346 


4.696 


1.499 


3.223 


3.516 


1.556 


4.886 


0.086 


0.217 


-0.20 


0.020 


0.048 


1.392 


5.509 


5.117 


1.604 


3.374 


4.078 


1.501 


5.125 


0.096 


0.234 


-0.10 


0.031 


0.060 


1.500 


5.659 


5.530 


1.708 


3.521 


4.646 


1.446 


5.355 


0.105 


0.252 


0.00 


0.043 


0.072 


1.607 


5.798 


5.943 


1.813 


3.666 


5.243 


1.392 


5.581 


0.116 


0.269 


0.10 


0.056 


0.086 


1.712 


5.924 


6.353 


1.917 


3.809 


5.875 


1.344 


5.805 


0.126 


0.287 


0.20 


0.071 


0.102 


1.816 


6.037 


6.767 


2.023 


3.953 


6.565 


1.302 


6.034 


0.136 


0.303 


0.30 


0.087 


0.119 


1.918 


6.132 


7.177 


2.128 


4.095 


7.308 


1.267 


6.266 


0.146 


0.320 


0.40 


0.108 


0.141 


2.018 


6.209 


7.590 


2.236 


4.241 


8.136 


1.238 


6.516 


0.155 


0.336 


0.50 


0.134 


0.169 


2.113 


6.260 


8.007 


2.344 


4.389 


9.057 


1.213 


6.779 


0.164 


0.352 


IFc/Hl 


Mg/; 


Fc527() 


Fc5335 


Fc54()6 


Fc57()9 


Fc57S2 


NaD 


H-,1 


l-l-iF 


H(),i 


H()/- 




-2.50 


0.411 


-0.137 


0.606 


0.227 


0.066 


-0.238 


1.001 


3.106 


2.968 


4.869 


3.830 




-2.40 


0.446 


-0.031 


0.592 


0.216 


0.089 


-0.201 


0.976 


2.966 


2.911 


4.807 


3.764 




-2.30 


0.487 


0.070 


0.584 


0.211 


0.113 


-0.163 


0.961 


2.881 


2.896 


4.782 


3.739 




-2.20 


0.531 


0.172 


0.582 


0.211 


0.137 


-0.126 


0.952 


2.837 


2.898 


4.779 


3.721 




-2.10 


0.578 


0.271 


0.585 


0.216 


0.163 


-0.088 


0.949 


2.838 


2.927 


4.801 


3.721 




-2.00 


0.634 


0.374 


0.597 


0.227 


0.191 


-0.050 


0.952 


2.784 


2.915 


4.771 


3.681 




-1.90 


0.696 


0.478 


0.617 


0.245 


0.220 


-0.010 


0.963 


2.729 


2.901 


4.738 


3.638 




-1.80 


0.771 


0.586 


0.648 


0.271 


0.252 


0.031 


0.984 


2.610 


2.847 


4.657 


3.565 




-1.70 


0.855 


0.696 


0.689 


0.305 


0.284 


0.073 


1.017 


2.481 


2.789 


4.564 


3.488 




-1.60 


0.951 


0.814 


0.742 


0.348 


0.317 


0.115 


1.065 


2.247 


2.670 


4.383 


3.356 




-1.50 


1.062 


0.938 


0.808 


0.399 


0.351 


0.159 


1.127 


1.909 


2.495 


4.119 


3.178 




-1.40 


1.196 


1.073 


0.890 


0.461 


0.388 


0.205 


1.208 


1.484 


2.271 


3.795 


2.969 




-1.30 


1.348 


1.208 


0.980 


0.528 


0.425 


0.250 


1.303 


1.046 


2.052 


3.479 


2.779 
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Table 5 — Continued 



Age = 12.0 



-1.20 


1.498 


1.338 


1.075 


0.597 


0.460 


0.293 


1.410 


0.637 


1.861 


3.160 


2.591 


-1.10 


1.678 


1.487 


1.189 


0.679 


0.498 


0.338 


1.533 


-0.051 


1.521 


2.619 


2.280 


-1.00 


1.878 


1.640 


1.314 


0.765 


0.538 


0.382 


1.665 


-0.868 


1.108 


1.979 


1.911 


-0.90 


2.113 


1.804 


1.456 


0.863 


0.581 


0.428 


1.815 


-1.796 


0.650 


1.296 


1.546 


-0.80 


2.352 


1.968 


1.603 


0.964 


0.625 


0.474 


1.976 


-2.687 


0.225 


0.650 


1.248 


-0.70 


2.597 


2.138 


1.764 


1.075 


0.675 


0.521 


2.153 


-3.548 


-0.167 


-0.022 


0.996 


-0.60 


2.825 


2.294 


1.918 


1.182 


0.722 


0.565 


2.332 


-4.178 


-0.423 


-0.563 


0.834 


-0.50 


3.042 


2.429 


2.047 


1.279 


0.772 


0.605 


2.465 


-4.757 


-0.631 


-1.035 


0.701 


-0.40 


3.251 


2.568 


2.183 


1.383 


0.821 


0.652 


2.607 


-5.301 


-0.815 


-1.513 


0.573 


-0.30 


3.458 


2.723 


2.358 


1.501 


0.867 


0.705 


2.828 


-5.747 


-0.965 


-1.984 


0.468 


-0.20 


3.658 


2.871 


2.534 


1.619 


0.919 


0.756 


3.057 


-6.181 


-1.109 


-2.441 


0.364 


-0.10 


3.857 


3.016 


2.712 


1.737 


0.974 


0.806 


3.296 


-6.598 


-1.251 


-2.886 


0.259 


0.00 


4.064 


3.158 


2.893 


1.857 


1.032 


0.855 


3.542 


-7.004 


-1.392 


-3.331 


0.152 


0.10 


4.272 


3.297 


3.078 


1.979 


1.089 


0.903 


3.804 


-7.392 


-1.536 


-3.771 


0.041 


0.20 


4.473 


3.437 


3.270 


2.105 


1.150 


0.952 


4.078 


-7.769 


-1.688 


-4.212 


-0.077 


0.30 


4.672 


3.575 


3.468 


2.234 


1.210 


0.999 


4.368 


-8.122 


-1.841 


-4.635 


-0.197 


0.40 


4.864 


3.716 


3.675 


2.367 


1.271 


1.051 


4.672 


-8.468 


-2.004 


-5.057 


-0.331 


0.50 


5.044 


3.864 


3.890 


2.500 


1.334 


1.105 


4.996 


-8.791 


-2.165 


-5.469 


-0.481 



Note. — The entire data of Table|5]are available at |http://web.yonsei.ac.kr/cosiiiic/dataA^EPS.htm 



Table 6. LICK absorption indices of YEPS simple stellar population model for [a/Fe]=03 (fitting 

functions of W94 and W97). 



Age = 12.0 



[Fe/H] 


CNi 


CN2 


Ca4227 


G4300 


Fe4383 


Ca4455 


Fe4531 


C24668 


H/J 


Fe5015 


Mgi 


Mg2 


-2.50 


-0.125 


-0.084 


0.235 


0.122 


0.611 


0.184 


0.371 


0.242 


2.561 


0.930 


0.008 


0.027 


-2.40 


-0.120 


-0.077 


0.274 


0.264 


0.570 


0.199 


0.500 


0.116 


2.576 


1.033 


0.008 


0.029 


-2.30 


-0.117 


-0.071 


0.298 


0.372 


0.529 


0.215 


0.618 


0.010 


2.601 


1.140 


0.008 


0.031 


-2.20 


-0.114 


-0.067 


0.322 


0.475 


0.502 


0.231 


0.733 


-0.080 


2.626 


1.251 


0.008 


0.033 


-2.10 


-0.110 


-0.062 


0.356 


0.606 


0.505 


0.252 


0.855 


-0.157 


2.634 


1.367 


0.009 


0.037 


-2.00 


-0.108 


-0.059 


0.399 


0.738 


0.517 


0.275 


0.977 


-0.214 


2.644 


1.484 


0.011 


0.042 


-1.90 


-0.106 


-0.056 


0.450 


0.872 


0.549 


0.301 


1.101 


-0.248 


2.654 


1.607 


0.014 


0.049 


-1.80 


-0.104 


-0.054 


0.507 


1.016 


0.602 


0.332 


1.230 


-0.265 


2.662 


1.738 


0.017 


0.057 


-1.70 


-0.103 


-0.052 


0.568 


1.164 


0.674 


0.366 


1.358 


-0.257 


2.666 


1.875 


0.021 


0.066 


-1.60 


-0.101 


-0.050 


0.640 


1.345 


0.777 


0.407 


1.493 


-0.222 


2.644 


2.030 


0.025 


0.076 


-1.50 


-0.100 


-0.049 


0.704 


1.482 


0.874 


0.447 


1.616 


-0.163 


2.646 


2.185 


0.029 


0.087 


-1.40 


-0.097 


-0.048 


0.775 


1.696 


1.023 


0.498 


1.754 


-0.095 


2.624 


2.362 


0.033 


0.099 


-1.30 


-0.090 


-0.044 


0.863 


2.032 


1.230 


0.564 


1.911 


0.019 


2.526 


2.569 


0.037 


0.112 


-1.20 


-0.079 


-0.038 


0.972 


2.481 


1.474 


0.647 


2.089 


0.179 


2.348 


2.805 


0.041 


0.126 


-1.10 


-0.067 


-0.031 


1.104 


2.994 


1.752 


0.744 


2.279 


0.400 


2.137 


3.068 


0.046 


0.142 


-1.00 


-0.059 


-0.026 


1.275 


3.464 


2.074 


0.847 


2.470 


0.659 


1.956 


3.348 


0.053 


0.158 


-0.90 


-0.055 


-0.025 


1.480 


3.836 


2.386 


0.950 


2.648 


0.922 


1.815 


3.622 


0.060 


0.174 


-0.80 


-0.053 


-0.025 


1.638 


4.089 


2.626 


1.036 


2.796 


1.209 


1.728 


3.864 


0.067 


0.190 


-0.70 


-0.050 


-0.023 


1.786 


4.327 


2.793 


1.117 


2.932 


1.489 


1.650 


4.129 


0.074 


0.206 


-0.60 


-0.041 


-0.014 


1.937 


4.583 


2.916 


1.200 


3.070 


1.796 


1.579 


4.448 


0.079 


0.222 


-0.50 


-0.034 


-0.006 


2.124 


4.766 


3.183 


1.291 


3.229 


2.240 


1.529 


4.705 


0.087 


0.242 


-0.40 


-0.028 


-0.001 


2.299 


4.942 


3.432 


1.375 


3.377 


2.642 


1.480 


4.940 


0.095 


0.262 


-0.30 


-0.022 


0.005 


2.466 


5.121 


3.687 


1.461 


3.526 


3.066 


1.436 


5.182 


0.104 


0.280 


-0.20 


-0.014 


0.012 


2.647 


5.315 


3.918 


1.544 


3.682 


3.513 


1.391 


5.443 


0.113 


0.298 


-0.10 


-0.006 


0.020 


2.815 


5.503 


4.157 


1.628 


3.838 


3.984 


1.348 


5.709 


0.121 


0.315 


0.00 


0.004 


0.030 


2.966 


5.684 


4.423 


1.717 


3.996 


4.488 


1.308 


5.986 


0.129 


0.331 


0.10 


0.014 


0.039 


3.079 


5.824 


4.721 


1.806 


4.144 


4.988 


1.264 


6.246 


0.136 


0.346 


0.20 


0.027 


0.051 


3.169 


5.951 


5.055 


1.902 


4.290 


5.529 


1.226 


6.519 


0.142 


0.359 


0.30 


0.042 


0.067 


3.241 


6.056 


5.425 


2.004 


4.434 


6.111 


1.192 


6.797 


0.147 


0.372 


0.40 


0.059 


0.084 


3.294 


6.131 


5.785 


2.102 


4.558 


6.760 


1.168 


7.062 


0.152 


0.383 


0.50 


0.079 


0.105 


3.349 


6.172 


6. 1 35 


2.199 


4.667 


7.486 


1.150 


7.301 


0.159 


0.396 


IFc/Hl 


Mg/; 


Fc527() 


Fc5335 


Fc5406 


Fc57()9 


Fc57S2 


NaD 


H-M 


l-l-iF 


H(),i 


Hdy 




-2.50 


0.556 


-0.050 


0.649 


0.237 


0.080 


-0.229 


1.032 


2.945 


2.886 


4.603 


3.685 




-2.40 


0.611 


0.060 


0.636 


0.228 


0.104 


-0.189 


1.033 


2.870 


2.870 


4.604 


3.656 




-2.30 


0.672 


0.165 


0.627 


0.222 


0.127 


-0.151 


1.034 


2.843 


2.887 


4.631 


3.655 




-2.20 


0.737 


0.269 


0.622 


0.221 


0.152 


-0.113 


1.042 


2.821 


2.903 


4.657 


3.652 




-2.10 


0.812 


0.379 


0.628 


0.228 


0.178 


-0.074 


1.060 


2.743 


2.879 


4.642 


3.614 




-2.00 


0.895 


0.489 


0.641 


0.242 


0.207 


-0.035 


1.093 


2.656 


2.848 


4.632 


3.576 




-1.90 


0.992 


0.602 


0.665 


0.265 


0.237 


0.006 


1.142 


2.553 


2.808 


4.616 


3.536 




-1.80 


1.101 


0.718 


0.699 


0.295 


0.268 


0.049 


1.207 


2.441 


2.760 


4.592 


3.482 




-1.70 


1.226 


0.834 


0.741 


0.333 


0.299 


0.091 


1.288 


2.296 


2.702 


4.543 


3.426 




-1.60 


1.376 


0.959 


0.797 


0.379 


0.332 


0.134 


1.386 


2.070 


2.598 


4.431 


3.337 




-1.50 


1.539 


1.075 


0.856 


0.428 


0.363 


0.176 


1.496 


1.921 


2.549 


4.371 


3.297 




-1.40 


1.709 


1.197 


0.924 


0.480 


0.394 


0.218 


1.612 


1.671 


2.442 


4.205 


3.172 




-1.30 


1.903 


1.331 


1.004 


0.540 


0.425 


0.259 


1.739 


1.134 


2.182 


3.777 


2.898 
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Table 6 — Continued 



Age= 12.0 



-1.20 


2.133 


1.479 


1.100 


0.608 


0.460 


0.302 


1.879 


0.326 


1.751 


3.130 


2.479 


-1.10 


2.419 


1.639 


1.210 


0.687 


0.499 


0.347 


2.036 


-0.636 


1.239 


2.418 


2.024 


-1.00 


2.744 


1.804 


1.330 


0.774 


0.544 


0.392 


2.219 


-1.593 


0.756 


1.797 


1.680 


-0.90 


3.069 


1.960 


1.449 


0.862 


0.592 


0.437 


2.416 


-2.361 


0.377 


1.302 


1.459 


-0.80 


3.364 


2.088 


1.554 


0.941 


0.636 


0.477 


2.618 


-2.823 


0.168 


0.980 


1.355 


-0.70 


3.650 


2.199 


1.639 


1.013 


0.682 


0.511 


2.793 


-3.220 


-0.005 


0.733 


1.287 


-0.60 


3.929 


2.301 


1.711 


1.084 


0.733 


0.550 


2.929 


-3.628 


-0.172 


0.510 


1.229 


-0.50 


4.200 


2.428 


1.832 


1.176 


0.778 


0.600 


3.175 


-3.889 


-0.284 


0.287 


1.194 


-0.40 


4.446 


2.546 


1.950 


1.261 


0.828 


0.646 


3.422 


-4.131 


-0.399 


0.085 


1.147 


-0.30 


4.676 


2.666 


2.074 


1.349 


0.882 


0.691 


3.683 


-4.381 


-0.513 


-0.138 


1.098 


-0.20 


4.893 


2.782 


2.191 


1.433 


0.939 


0.736 


3.970 


-4.593 


-0.633 


-0.312 


1.057 


-0.10 


5.114 


2.899 


2.318 


1.519 


0.997 


0.780 


4.276 


-4.804 


-0.759 


-0.498 


1.012 


0.00 


5.333 


3.024 


2.457 


1.614 


1.057 


0.825 


4.593 


-5.038 


-0.899 


-0.722 


0.954 


0.10 


5.547 


3.150 


2.612 


1.716 


1.113 


0.871 


4.915 


-5.280 


-1.045 


-0.981 


0.882 


0.20 


5.747 


3.284 


2.781 


1.826 


1.171 


0.919 


5.239 


-5.561 


-1.203 


-1.304 


0.794 


0.30 


5.929 


3.426 


2.965 


1.946 


1.231 


0.971 


5.566 


-5.881 


-1.376 


-1.691 


0.685 


0.40 


6.089 


3.563 


3.147 


2.064 


1.285 


1.019 


5.876 


-6.171 


-1.536 


-2.071 


0.571 


0.50 


6.214 


3.698 


3.326 


2.179 


1.336 


1.067 


6.168 


-6.443 


-1.694 


-2.448 


0.439 



Note. — The entire data of Table[6]are available at |http://web.yonsei.ac.kr/cosmic/dataA^EPS.htni| 
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Table 7. LICK absorption indices of YEPS simple stellar population model for [a/Fe]=0.6 (fitting 

functions of W94 and W97). 



Age = 12.0 



[Fe/H] 


CNi 


CN2 


Ca4227 


G4300 


Fe4383 


Ca4455 


Fe4531 


C24668 


H/J 


Fe5015 


Mgi 


Mg2 


-2.50 


-0.126 


-0.082 


0.207 


-0.015 


0.733 


0.206 


0.464 


0.409 


2.587 


1.073 


0.010 


0.031 


-2.40 


-0.123 


-0.076 


0.264 


0.092 


0.656 


0.218 


0.593 


0.288 


2.622 


1.179 


0.011 


0.035 


-2.30 


-0.119 


-0.071 


0.319 


0.216 


0.583 


0.227 


0.721 


0.156 


2.630 


1.284 


0.010 


0.037 


-2.20 


-0.115 


-0.066 


0.378 


0.340 


0.519 


0.239 


0.848 


0.043 


2.635 


1.393 


0.011 


0.041 


-2.10 


-0.113 


-0.062 


0.443 


0.460 


0.477 


0.254 


0.978 


-0.040 


2.651 


1.507 


0.013 


0.047 


-2.00 


-0.110 


-0.059 


0.518 


0.600 


0.473 


0.275 


1.114 


-0.092 


2.663 


1.635 


0.016 


0.055 


-1.90 


-0.107 


-0.055 


0.603 


0.769 


0.516 


0.306 


1.259 


-0.102 


2.657 


1.780 


0.020 


0.065 


-1.80 


-0.104 


-0.052 


0.696 


0.951 


0.589 


0.342 


1.408 


-0.083 


2.648 


1.938 


0.025 


0.077 


-1.70 


-0.102 


-0.049 


0.787 


1.119 


0.678 


0.381 


1.553 


-0.045 


2.657 


2.103 


0.030 


0.090 


-1.60 


-0.096 


-0.045 


0.896 


1.396 


0.827 


0.433 


1.716 


0.009 


2.605 


2.294 


0.035 


0.104 


-1.50 


-0.087 


-0.040 


1.025 


1.782 


1.015 


0.497 


1.896 


0.077 


2.481 


2.504 


0.040 


0.119 


-1.40 


-0.077 


-0.034 


1.175 


2.247 


1.215 


0.574 


2.089 


0.182 


2.295 


2.740 


0.045 


0.134 


-1.30 


-0.066 


-0.027 


1.363 


2.730 


1.433 


0.661 


2.293 


0.334 


2.094 


3.004 


0.051 


0.152 


-1.20 


-0.060 


-0.024 


1.631 


3.099 


1.637 


0.750 


2.501 


0.484 


1.926 


3.297 


0.057 


0.169 


-1.10 


-0.060 


-0.026 


1.892 


3.378 


1.791 


0.831 


2.677 


0.629 


1.802 


3.551 


0.064 


0.187 


-1.00 


-0.061 


-0.030 


2.120 


3.587 


1.880 


0.898 


2.821 


0.779 


1.714 


3.773 


0.071 


0.205 


-0.90 


-0.062 


-0.032 


2.347 


3.814 


1.887 


0.961 


2.956 


0.913 


1.621 


4.019 


0.077 


0.223 


-0.80 


-0.057 


-0.028 


2.585 


4.059 


1.863 


1.029 


3.099 


1.106 


1.533 


4.316 


0.082 


0.243 


-0.70 


-0.056 


-0.027 


2.838 


4.223 


1.981 


1.102 


3.251 


1.397 


1.473 


4.548 


0.089 


0.264 


-0.60 


-0.056 


-0.028 


3.088 


4.383 


2.098 


1.175 


3.400 


1.672 


1.413 


4.770 


0.099 


0.287 


-0.50 


-0.056 


-0.029 


3.334 


4.545 


2.225 


1.249 


3.548 


1.960 


1.357 


4.992 


0.108 


0.309 


-0.40 


-0.056 


-0.029 


3.602 


4.717 


2.329 


1.321 


3.702 


2.256 


1.301 


5.226 


0.119 


0.331 


-0.30 


-0.053 


-0.027 


3.857 


4.923 


2.452 


1.398 


3.866 


2.594 


1.257 


5.497 


0.128 


0.352 


-0.20 


-0.048 


-0.023 


4.093 


5.150 


2.583 


1.475 


4.037 


2.977 


1.225 


5.803 


0.137 


0.370 


-0.10 


-0.042 


-0.018 


4.295 


5.374 


2.748 


1.557 


4.209 


3.386 


1.202 


6.124 


0.144 


0.386 


0.00 


-0.033 


-0.011 


4.463 


5.599 


2.964 


1.644 


4.380 


3.808 


1.186 


6.458 


0.148 


0.399 


0.10 


-0.026 


-0.005 


4.561 


5.749 


3.237 


1.732 


4.534 


4.167 


1.155 


6.751 


0.151 


0.411 


0.20 


-0.016 


0.003 


4.645 


5.880 


3.567 


1.835 


4.696 


4.536 


1.121 


7.067 


0.152 


0.421 


0.30 


-0.003 


0.015 


4.746 


5.974 


3.971 


1.961 


4.874 


4.888 


1.079 


7.392 


0.156 


0.434 


[Fe/H] 


Mgh 


Fe5270 


Fe5335 


Fe54{)6 


Fe5709 


Fe5782 


NaD 




H7/r 


H5a 


H<5,r 




-2.50 


0.689 


0.031 


0.715 


0.261 


0.097 


-0.214 


1.054 


3.152 


3.079 


4.586 


3.776 




-2.40 


0.767 


0.139 


0.697 


0.250 


0.120 


-0.175 


1.078 


3.138 


3.103 


4.658 


3.783 




-2.30 


0.855 


0.239 


0.673 


0.236 


0.142 


-0.140 


1.097 


3.142 


3.098 


4.759 


3.772 




-2.20 


0.952 


0.339 


0.655 


0.228 


0.164 


-0.104 


1.127 


3.146 


3.086 


4.872 


3.764 




-2.10 


1.057 


0.444 


0.650 


0.231 


0.187 


-0.067 


1.176 


3.139 


3.081 


4.974 


3.756 




-2.00 


1.174 


0.559 


0.663 


0.246 


0.213 


-0.027 


1.250 


3.050 


3.048 


5.011 


3.722 




-1.90 


1.313 


0.687 


0.694 


0.275 


0.241 


0.016 


1.350 


2.853 


2.973 


4.950 


3.652 




-1.80 


1.471 


0.821 


0.739 


0.314 


0.271 


0.060 


1.471 


2.600 


2.883 


4.834 


3.567 




-1.70 


1.643 


0.952 


0.792 


0.357 


0.300 


0.104 


1.607 


2.381 


2.827 


4.732 


3.498 




-1.60 


1.838 


1.095 


0.855 


0.406 


0.330 


0.149 


1.755 


1.956 


2.653 


4.444 


3.306 




-1.50 


2.059 


1.243 


0.925 


0.458 


0.361 


0.192 


1.909 


1.296 


2.318 


3.947 


2.968 




-1.40 


2.324 


1.399 


1.004 


0.516 


0.395 


0.235 


2.072 


0.501 


1.891 


3.361 


2.559 




-1.30 


2.653 


1.562 


1.091 


0.581 


0.433 


0.280 


2.260 


-0.373 


1.428 


2.795 


2.166 




-1.20 


3.022 


1.725 


1.175 


0.649 


0.474 


0.324 


2.476 


-1.120 


1.041 


2.411 


1.910 




-1.10 


3.394 


1.857 


1.245 


0.709 


0.515 


0.363 


2.701 


-1.599 


0.785 


2.199 


1.784 
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Table 7 — Continued 



Age = 12.0 



-1.00 


3.751 


1.957 


1.296 


0.757 


0.553 


0.397 


2.925 


-1.834 


0.641 


2.155 


1 


.756 


-0.90 


4.110 


2.033 


1.322 


0.795 


0.594 


0.424 


3.117 


-2.035 


0.508 


2.199 


1 


.763 


-0.80 


4.468 


2.102 


1.340 


0.836 


0.638 


0.455 


3.278 


-2.223 


0.384 


2.280 


1 


.794 


-0.70 


4.802 


2.194 


1.405 


0.899 


0.680 


0.499 


3.528 


-2.273 


0.308 


2.318 


1 


.822 


-0.60 


5.116 


2.284 


1.476 


0.962 


0.729 


0.542 


3.790 


-2.324 


0.231 


2.357 


1 


.842 


-0.50 


5.416 


2.377 


1.556 


1.029 


0.782 


0.586 


4.061 


-2.400 


0.136 


2.376 


1 


.846 


-0.40 


5.712 


2.468 


1.635 


1.094 


0.838 


0.629 


4.362 


-2.445 


0.041 


2.443 


1 


.858 


-0.30 


5.978 


2.568 


1.722 


1.165 


0.901 


0.674 


4.680 


-2.539 


-0.072 


2.441 


1 


.853 


-0.20 


6.224 


2.673 


1.816 


1.237 


0.966 


0.719 


5.031 


-2.645 


-0.193 


2.396 


1 


.841 


-0.10 


6.460 


2.788 


1.925 


1.315 


1.031 


0.765 


5.392 


-2.789 


-0.319 


2.289 


1 


.816 


0.00 


6.677 


2.913 


2.049 


1.400 


1.098 


0.812 


5.760 


-2.995 


-0.462 


2.088 


1 


.767 


0.10 


6.894 


3.041 


2.198 


1.494 


1.154 


0.857 


6.127 


-3.209 


-0.605 


1.843 


1 


.701 


0.20 


7.111 


3.190 


2.378 


1.610 


1.217 


0.911 


6.505 


-3.488 


-0.767 


1.523 


1 


.617 


0.30 


7.330 


3.364 


2.599 


1.756 


1.292 


0.980 


6.918 


-3.864 


-0.960 


1.088 


1 


.499 



Note. — The entire data of Table|2]are available at |http://web.yonsei.ac.kr/cosmic/data/YEPS .htm | 



